Finite-Size Effects on Critical Diffusion and Relaxation Towards Metastable 

Equilibrium 
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We present the first analytic study of finite-size effects on critical diffusion above and below T c 
of three-dimensional Ising-like systems whose order parameter is coupled to a conserved density. 
We also calculate the finite-size relaxation time that governs the critical order-parameter relaxation 
towards a metastable equilibrium state below T c . Two new universal dynamic amplitude ratios at 
T c are predicted and quantitative predictions of dynamic finite-size scaling functions are given that 
can be tested by Monte-Carlo simulations. 



PACS numbers: 64.60.Ht, 75.40.Gb, 75.40.Mg 

The dissipative critical dynamics of bulk systems with 
a non-conserved order parameter are fairly well under- 
stood. Depending on whether the order parameter is 
governed by purely relaxational dynamics or whether it 
is coupled to a hydrodynamic (conserved) density such 
systems belong to the universality classes of models A 
or C . The fundamental dynamic quantities of these 
systems are the relaxation and diffusion times which di- 
verge as the critical temperature T c is approached. 

For finite systems, these times are expected to become 
smooth and finite throughout the critical region and to 
depend sensitively on the geometry and boundary condi- 
tions. These finite-size effects are particularly large in 
Monte Carlo (MC) simulations of small systems. On 
a qualitative level, they can be interpreted on the ba- 
sis of phcnomenological finite-size scaling assumptions. 
For a more stringent analysis the knowledge of the shape 
of universal finite-size scaling functions is necessary. So 
far there exist reliable theoretical predictions on finite- 
size dynamics in three dimensions only on two relaxation 
times T\ and determining the long-time behavior of the 
order parameter and the square of the order parameter 
HU . No analytic work exists, to the best of our knowl- 
edge, on the important universality class JjJ of diffusive 
finite-size behavior near T c . This is of relevance, e.g., to 
magnetic systems with mobile impurities JjJ , to binary 
alloys with an order-disorder transition to uniaxial 
antiferromagnets or to systems whose order parame- 
ter is coupled to the conserved energy density jj],§. 

In this Letter we present the first prediction of the 
finite-size scaling function for the critical diffusion time 
of three-dimensional systems above and below T c . Fur- 
thermore we shall present the analytic identification and 
quantitative calculation of a new leading relaxation time 
that governs the critical order-parameter relaxation to- 
wards a metastable equilibrium state of finite systems be- 
low T c . Our predictions contain no adjustable parameters 
other than two amplitudes of the bulk system. 

We start from model C 0], i.e., from the relax- 
ational and diffusive Langevin equations for the one- 



component order-parameter field <p(x, t) and for the den- 
sity p(x, t) = (p) + m(x, t) in a finite volume V, 



dy(x, t) 

dt 
<9m(x, t) 



= -r„ 



SH 



= An V 2 



SH 



+ e„(x,t), 



+ 9 m (x,t) , 



dt 5m(x, t) 

H = Iv ddx [ \ T ^ 2 + U^<P) 2 + <W 4 
+ \m 2 + 7oto<^ 2 - h m] 



(1) 
(2) 

(3) 



where Q v and O m are Gaussian (S-correlated random 
forces. We consider an equilibrium ensemble near T c (p) 
at fixed p = V^ 1 J v d d x p(x, t). This corresponds to the 
experimental situation of keeping the conserved quan- 
tity (e.g., number of impurities) fixed when changing 
the reduced temperature t = [T - T c (p)]/T c (p). The 
latter enters through to. Because of (p) = p we have 
(m) =m = 0. Eqs. (0) - (||) describe the dynamics of re- 
laxational and diffusive modes that are coupled through 
7o- We are interested in the long-time behavior of the 
diffusive modes above, at and below T c as well as in the 
order-parameter relaxation on an intermediate time scale 
below T c . We shall begin with the diffusive modes. For 
simplicity and for the purpose of a comparison with MC 
simulations, we assume cubic geometry, V = L d , with 
periodic boundary conditions. 

For the bulk system, the diffusion constants D^it) 
above and below T c appear in the small k limit of the 
long-time behavior of the correlation function 

C n (k,t,t) = l/-> k (i)7i_ k (0)) - exp [-D ± (t)k 2 t] (4) 

where n\^{t) — nik(t) + c n (k)?pk(t) is an appropriate lin- 
ear combination of mk(i) = j v d d x m(x, i)e~ jkx and 
V>k(*) = J v d d x [<p(x, t) - {ip)\ e- 2kx with c n (k) = ~c n k 2 + 
0(k 4 ). The coefficient c„ can be identified by linearizing 
Eqs. (@) - ©. Above T c , c„ = because of (ip) = 0. 
At T c , the long-time behavior of C n is non-exponential 
(power law) for the bulk system. 

For the finite system, the coefficient c n (k) is modified 
[via the replacement (ip) — > M as defined in Eqs. (0) 
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and (||) below] and the long-time behavior of C n remains 
exponential, 

C n (k,t,L,t) ~ exp[-Q„(M,L)t] , (5) 

even in the non-hydrodynamic region at bulk T c where 
the small- k approximation is no longer justified. As a 
conceptual complication there exists a smallest nonzero 
value fc^ in = 4ir 2 /L 2 of k 2 which prevents us to perform 
the limit k — > for the finite system. Therefore we need 
to derive the finite-size scaling function for Cl n (k,t, L) at 
finite k. Nevertheless we may define an effective diffu- 
sion time td = ^ n (27rL _1 , t, L)^ 1 or a diffusion constant 
D = il n /k 2 at k — k rn i n — 2ir/L of the finite system by 

D(t,L) = (2n)- 2 L 2 r2„(2 7 rL~ 1 ,f,L) (6) 

which interpolates smoothly between the bulk result 
D(t, oo) = D ± {t) above and below T c (Fig. Q). 
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FIG. 1. Diffusion constant D(t,L)/A+ for L = 80£o_(solid 
line) vs t, and of the scaling function Jd(x), Eq. (113), vs 
x = tL [1 ' a)/v (solid line), with L in units of f - The dashed 
lines represent the bulk diffusion constants D ± (t)/A~] D . 

In the spirit of finite-size theory |3|J^] we decompose 
</j(x, t) = Mo + i5(/?(x, t) with the zero-mode average 

/oo />oo 
dMM 2 e- Ho / / dMe- Ho (7) 
-oo J —oo 

where ff (M) = L d (±T M 2 + u M 4 ) is the = part of 
H, with M = V~ l J v d d x (p. For the finite system, the 
quantity Mq(t , L) is non-zero for all T and interpolates 
smoothly between T > T c and T < T c . Linearization of 
Eqs. (|f|) - (||) with respect to S<pk(t) and rok(i) leads to 

cn(k) = (wr 1 ^ - [(6 -) 2 + ^o7o 2 fc 2 ] * } , (8) 

Q n (k,t,L) = ±A {6+ - [(6 -) 2 +w^k 2 Y } , (9) 
6±(fc) = w (to + l2u a M 2 + k 2 ) ± fc 2 , (10) 



with iuo = To/Ao and 70 = 470-M0. 

An application of these unrenormalized expressions 
to the critical region requires us to turn to the renor- 
malized theory. The strategy of the field-theoretic 
renormalization-group (RG) approach at d — 3 dimen- 
sions is well established in bulk statics || and dynam- 
ics H and has been successfully applied recently to the 
model-A finite-size dynamics || . The details of its appli- 
cation to model C will be given elsewhere [^0|. Here we 
only present the asymptotic finite-size scaling form 

n n (k,t,L) = L-* UiL^'v ,kL) (II) 

as derived from Eqs. (|) and (|) - ([k]), with the dynamic 
critical exponent z = 2 + ajv jl. The scaling function 
reads in three dimensions 

f n (x, K ) = A n l a / v {b + - [b 2 _ + w*c*t l ' 2 K H 2 {y)] k } , 
b±(x,K) = w*[£(x) 2 + k 2 ]±k 2 , 
l{xf' 2 = ^u*f/ 2 {y{x) + 120 a (jf(a:))} , 
y(x) = (ATTu*)- 1/2 £{x) {z/2) - {1 - a)/w x , 
MV) = Uo°°ds s 2 e-^ 2 ^)/{^ds e-^^l 

with c* = 16(7*) 2 (4 7 r/{t*) 1 /2 and x = x ^ {1 ~ a)/v . 
This yields the scaling form D(t,L) = L 2 ^ z fn(x) for 
the diffusion constant, Eq. (||), with 

f D (x) = (2 7 r)- 2 /„( 2 ;,27r). (12) 

The static parameters are Jl u* = u* + (7*) 2 /2 and 
( 7 *) 2 = a[4i/B(u*)]~ 1 with |L1| u* = 0.0404 and B(u*) = 
0.502 in three dimensions. For v and a we take 0.6335 
and 0.100 p2| . The dynamic parameter is w* = 1 in 
one-loop order. The two non-universal bulk amplitudes 
£o(p) and A n = ^A^q^ 2 are defined by the asymptotic 
behavior £ = foi -1 ^ 1-00 and D+{t) = A+t^- 2 ^^ 1 -^ 
of the correlation length and diffusion constant at fixed 
p above T c . The exponent vj{\ — a) instead of v is due 
to Fisher renormalization fll3|| . 

The solid line in Fig. [l] shows D(t, L)jA + D vs t for the 
example L = 80£o- The same line represents /d{x) vs x 
(top scale). For comparison the bulk limits (dashed 
lines) are also shown, with A D /A^ = 2 Q /( 1 ~ Q )(1 + 
\"f* 2 /u*)^ 1 = 0.55. We expect the accuracy of these 
results to be of O(10%). These predictions can be tested 
by MC simulations, after adjusting £o(p) and in the 
bulk region i»l above T c . 

In addition to the finite-size effect on the diffusive 
modes there exists an interesting finite-size effect on the 
relaxational modes below T c that has so far not been in- 
vestigated analytically. It is well known that no sponta- 
neous symmetry breaking can take place in finite systems 
below T c because of ergodicity. For Ising-like systems, er- 
godicity implies a " tunneling" between metastable states 
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of opposite orientation of the magnetization as observed 
in MC simulations ]T^-[T^| . On an intermediate time scale 
t < t x (V), however, the magnetization does not change 
sign and its magnitude relaxes towards a finite value that 
characterizes such a metastable state |E3|. This relax- 
ation process is important for large systems since the 
crossover time t x (L) is expected to grow with the size L 
as ~ L z where z is the dynamic critical exponent. This 
process occurs both in model A and model C, therefore 
we confine ourselves to the simpler model A in the follow- 
ing. We stress that the relaxation process for t < t x {L) is 
fundamentally different from the ultimate long-time be- 
havior for t ^> t x (L) studied previously HQ]. 

Model A is defined by Eq. (|l|) where H is replaced by 

H v = f d d x [ \ r cp 2 + I (Vtp) 2 + u cp 4 ] . (13) 
Jv 

We consider the time-dependent spatial average 
M(t,L) = L~ d J v d d xip(x,t) . We are primarily inter- 
ested in the long-time behavior of the equilibrium cor- 
relation function (M(t, L)M(0, L)) = C(t,L) for d = 3. 
For the bulk system, this behavior is 

C(t,oo)~^+exp(-t/r+), (14) 
C(t, oo) - M% ~ A 6 " exp(-t/r 6 -) (15) 

above and below T c where M sp — lim t ^ 00 (M(t, oo)} is 
the spontaneous order parameter and are the bulk re- 
laxation times. For the finite system above T c , the leading 
time dependence is still a single exponential ~ c 1 e _ *^ ri ' i - ) 
with a relaxation time T\(L) whose finite-size scaling 
function is known both analytically ||]4| and numerically 
||[l6|. In particular, lirriL-too t\{L) — rjf~. 

For the finite system below T c , however, the situation 
is more complicated and considerably less well explored. 
MC simulations JlJ] and phenomenological considera- 
tions suggest that there should exist an L dependent 
generalization t~(L) of t 6 ~~ , with lim£^ 00 r _ (i) = T b ~ , 
which should describe (i) the exponential relaxation of 
{M(t, L)) towards a metastable finite value on an inter- 
mediate time scale t < t x (L) before tunneling sets in, 
and (ii) a corresponding exponential decay of C(t,L) on 
this time scale. The question arises whether and how this 
important relaxation time t~{L) can be identified ana- 
lytically within models A and C. This question was left 
unanswered in the previous literature. In particular, nei- 
ther ti(L) nor T2(L) as calculated previously ||, can be 
identified with t~(L). [Below T c ,ti describes the decay 
of (M(t,L)) and of C(t,L) towards zero for t 3> t x (L) 
due to tunneling processes, and ti describes the decay 
of (M(t) 2 ) and of {M(t) 2 M(0) 2 ) towards {M 2 ) eq and 
{M 2 ) 2 q , respectively, for t ^> t x (L).] In the following we 
establish an analytic identification of t~ (L) and present 
the first quantitative prediction for its finite-size scaling 
behavior. 

To elucidate the main features we first neglect the 



inhomogeneous fluctuations cr(x, t) — tp(x.,t) — M(t). 
Then Eq. (|l|) is equivalent to the Fokker-Planck equation 
dP{M,t)/dt = -C P(M,t) for the probability distribu- 
tion P(M, t) with the operator 

T 9 f dH (M) d s 

c "-'T^dMy~^M~ + dM) (16) 

where H (M) = L d (\r Q M 2 + u M 4 ) . It is well known 
]l7t that C(t, L) is determined by the eigenvalues tk and 
eigenfunctions <pk(M) of Co according to 

oo 

C(t,L) = J2ck{L)eM-t/T k {L)], t>0 (17) 
k=i 

with cfc(L) = [JZo dM M MM)] 2 and r k (L) = e~\ 
E{) = < ei < 62 . . . By symmetry, c k = for even 
values of k. Below T c , T\(L) diverges in the bulk limit 
and lim£_ +00 cie _t / Tl = M 2 p becomes time-independent, 
thus an analysis of the k = 3 term in Eq. (]l7]) becomes 
indispensable. From the spectrum of Cq |l7j we find a de- 
generacy for k = 3 and k — 5 in the bulk limit for ro < 0. 
This requires to take the k = 5 term into account as well. 
We have found, however, that the coefficient c 5 vanishes 
in the bulk limit below T c whereas C3 remains finite. For 
finite L near T c , T5 is well separated from T3, as shown 
below. Thus it suffices to describe the time dependence 
of C(t, V) on intermediate time scales t ~ 0(t~(L)) and 
0(t~(L)) < t < 0(ti(L)) as well as on the long-time 
scale t ^ ti (L) as 

C(t,L) ~ci(L)e- t/Tl(L) + c 3 (L)e- t/T3(L) (18) 

where ci(oo) = M 2 p and €3(00) = below T c . In 
particular we arrive at the desired identification 

t-(L) = t 3 (L), lim r 3 (L) = t- . (19) 

L — >oo 

We conclude that, although t^{L) represents only a sub- 
leading relaxation time above T c , t^{L) governs the lead- 
ing time dependence ofC(t,L) of large finite systems be- 
;i :g.|i. 

These results yield also the key to the interpretation 
of t^(L) as the relaxation time governing the approach 
of the non- equilibrium quantity (M (t, L)) towards a 
metastable finite value before M(t, L) starts to change 
sign. This interpretation is based on the fact (lj] that the 
leading relaxation times of (M(t,L)) are determined by 
the same eigenvalues of Cq as the long-time behavior of 
the equilibrium correlation function C, i.e., 

(M(t, L)) - ci(L) e" t/Tl(L) + 5 3 (i) e~* /r3(i) . (20) 

The basic difference between Eqs. ( p0| ) and ( |l8| ) is that 
the coefficients Ck depend on the initial (non-equilibrium) 
state. 
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We proceed by presenting the results of a quantita- 
tive calculation of t 3 (L) and t 5 (L) including the effect 
of the inhomogeneous fluctuations ct(x) to one-loop or- 
der. This calculation is parallel to that performed pre- 
viously || and is expected to be as reliable as the pre- 
vious results pi . It is based on the Fokker-Planck equa- 
tion dP(M,t)/dt = -dP{M,t) where d has the same 
structure as Co, Eq. (16), but with tq, uq, Tq replaced by 
(positive) effective parameters j-J) , uj^ , T^f* ||f7],|l9| . 
In terms of the eigenvalues /^(k) and ^(k) of the equiv- 
alent Schrodinger equation jlTj we determine the relax- 
ation times T3 and T5 as 



= ( 2 r^)- 1 i rf / 2 (u^)- 1/2 /*(*) 



(21) 



with k = ^Tq L d / 2 (u^/f )~ 1 / 2 . In the asymptotic region 
the field-theoretic RG approach at d — 3 yields 
the finite-size scaling form Tj = L z fi(x),i = 3,5, with 
the scaling variable x = tL 1 ^, t = (T - T c )/T c . The an- 
alytic expressions for fi{x) are analogous to those given 
previously (|] and will be given elsewhere E3j. At T c we 
predict the universal ratios ti/ts = 8.5 and T3/T5 = 2.3. 
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FIG. 2. Relaxation times n(t,L)/A+ b for L = 80a (solid 
lines) vs i= (T— T c )/T c , and of their scaling functions fi(x) 
vs x — t L 1 ^ (solid lines), with L in units of the lattice con- 
stant a. Dashed lines: bulk relaxation times T^(t)/A^ b . 

The results are shown in Fig. |[ For an application 
to the Ising model we have taken £ /a = 0.495 fl| 
where a is the lattice spacing. The relaxation times Tj 
in Fig. are normalized to the bulk amplitude A^ b of 
t+ = A+ b t- yz 7 z = 2.04 (dashed line above T c ). Be- 
low T c , our theory yields the expected (2^^] exponen- 
tial decay, Eq. (|l5|), for the d — 3 bulk system, in dis- 
agreement with Rcf. p2]. The dashed line below T c rep- 



resents the bulk relaxation time r b = A Tb \t\ vz with 



A~ b /A+ = 2— (1 + fO/(l + 
mensions. Unlike for T\ and T2 
presently available for T3. 
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0.26 in three di- 
, no MC data are 



In summary we have presented the first quantitative 
predictions for the finite-size effects on critical diffusion 
and order-parameter relaxation towards metastable equi- 
librium in three-dimensional systems near T c . It would 
be interesting to test the predicted universal ratios ti/t%, 
T3/T5 and the finite-size scaling functions /d(x) and fi(x) 
(Figs. 1 and 2) by MC simulations. This appears to be 
within reach of present simulation techniques ]23] ] . 

Support by SFB 341 der Deutschen Forschungsgemein- 
schaft and by NASA is acknowledged. 
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